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Abstract 

The mass shift Am p and width broadening Ar p of /9-mesons produced in heavy 
ion collisions is estimated using general formule which relate the in-medium mass 
shift of a particle to the real part of the forward scattering amplitude Hef(E) of this 
particle on constituents of the medium and Ar to the corresponding cross section. 
It is found that the mass increases by some tens of MeV but, more importantly, the 
width becomes large, increasing by several hundred MeV at beam energies of a few 
GeVA and by twice that amount at beam energies of about a hundred GeVA. 
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1 Introduction. 



The problem of how the properties of hadrons change in hadronic or nuclear matter 
in comparison to their free values has attracted a lot of attention recently. Among these 
properties of immediate interest are the in-medium particles mass shifts and width broad- 
enings. Different models, as well as model independent approaches, were used to calculate 
these effects both at finite temperature and finite density. For a review see It is clear 
on physical grounds that the in-medium mass shift and width broadening of a particle 
is only due to its interaction with the constituents of the medium. Thus one can use 
phenomenological information on this interaction to calculate the mass shifts. In a recent 
paper 0, two of us have argued that the mass shift of a particle in medium can be re- 
lated to the forward scattering amplitude f(E) of this particle on the constituents of the 
medium. Written in the medium rest frame it is: 

Am(E) = —2-K—Ref(E) . (1) 
m 

Here m is the vacuum mass of the particle, E is its energy in the rest frame of the con- 
stituent particle, and p is the density of consituents. The normalization of the amplitude 
corresponds to the standard form of the optical theorem, 

ka = Anlmf(E) , (2) 
where k is the particle momentum. The width broadening is given by 

AT(E) = — ka(E) . (3) 
m 

The domain of applicability of Eqs. (1) and (3) was discussed in [f|. Briefly: 

• The particle's wavelength A must be much less than the mean distances between medium 
constituents d: A = k^ 1 <C d. This means that the particle's momentum k must be larger 
than a few hundred MeV. 

• The particle's formation length If ~ (E/m)/m c } iai -, with m c har ~ mp, must be less than 
the nucleus radius R. 

• Ref(E), which enters Eq. (1), must satisfy the inequality | Re/ |< d. 

• The main part of the scattering proceeds through small angles, < 1. Only in this 
case is the optical analogy on which Eqs. (1) and (3) are based correct. 

Equations (1) and (3) are correct also when the medium constituents have some mo- 
mentum distributions, such as Fermi-Dirac or Bose-Einstein distributions for finite tem- 
peratures and chemical potentials. In these cases averaging over the constituents' mo- 
mentum distributions must be performed on the right sides of Eqs. (1) and (3). These 
equations were derived in [2] based on simple quantum mechanical arguments and the op- 
tical analogy. This approach allows one to formulate in an explicit way the applicability 
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conditions presented above. When the medium is a gas in thermal equilibrium the equiv- 
alent to Eqs. (1) and (3) can be derived on the basis of thermal field theory. References 
]3|, [| give a few examples and reference H gives a relativistic field-theoretic derivation. 

In most of the papers on the in-medium hadron mass shifts the hadrons were considered 
at rest. As seen from Eq. (1) this restriction is not necessary theoretically. It is desirable 
to have theoretical predictions in a broad energy range since it extends the possibilities of 
experimental investigation. As discussed in |J for the cases of p or 7r-mesons embedded 
in nuclear matter the energy dependence of the mass shifts is rather significant at low 
energies where resonances dominate. 

We estimate the mass shift and width broadening in the case of p-mesons produced 
in heavy ion collisions. The most interesting case is that of the p° which can be observed 
through the decay p° — > e + e~ or p + p~. We will assume that p- mesons are formed in the 
last stage of the evolution of hadronic matter created in course of a heavy ion collision 
when the matter can be considered as an almost noninteracting gas of pions and nucleons. 
(We will neglect the admixture of kaons and hyperons, which is known to be small 0, as 
well as heavy resonances.) This stage occurs late in the collision when the total density 
of nucleons and pions is of the order of the normal nucleon density in a nucleus. The 
description of nuclear matter as a noninteracting gas of nucleons and pions cannot be 
considered as a very good one so it is clear from the beginning that our results may be 
only semiquantitave. The main ingredients of our calculation are pit and pN forward 
scattering amplitudes and total cross sections as well as the values of nucleon and pion 
densities. 

In this paper we consider central heavy ion collisions and assume that nucleon and 
pion momentum distributions in the gas are just the momentum distributions measured 
experimentally in such collisions. The case where nucleons and pions are assumed to exist 
in a state of equilibrium at fixed temperature and chemical potential will be considered 
in a subsequent publication [8]. 

2 Calculation of pN and pn cross sections and forward 
scattering amplitudes. 

Let us first focus on the amplitudes and cross sections. To determine these quantities we 
use the following procedure. At low energies we saturate the cross sections and forward 
scattering amplitudes by resonance contributions. At high energies we determine a p ^ and 
a p7T from cr 7 7v and <r 77r using the vector dominance model (VDM). The cross section ct 7 tv 
is well known experimentally [|§, Re/ 7 Ar is determined from the dispersion relation, and 
a 77r and Re/ 77r can be found by the Regge approach. Since VDM allows one to find only 
the cross sections of transversally polarized p-mesons we restrict ourselves to this case. 
As was shown in [2] , when E p > 2 GeV, Am and Ar for longitudinal p-mesons are much 
smaller than for transversal ones in nuclear matter. At zero p-meson energy, Am and Ar 
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for transverse and longitudinal p-mesons are evidently equal. In the case of scattering on 
a low temperature pion gas they are comparable [10]. Therefore our results should be 
multiplied by a factor ranging between 1 and 1.5 for unpolarized p-mesons. 

To estimate Re/ p7r (s) at low energy we write, in the center of mass (cm.) frame: 

Ref ( S ) = -Y FF .J: B R T R (^-m R ) 

where y/s is the total cm. energy, m R and T R are the mass and total width of the 
resonance, B R is the branching ratio of its decay into np and q cm is the center of mass 
momentum 



Qcm = \J[s - (m p + m n ) 2 ][s - (m p - m 7T ) 2 ]/2y/s . (5) 

Here F s is the spin factor and is the isospin factor. The latter is equal to (1/2) x (2/3) = 
1/3 for I R = 1. The first factor reflects the fact that we are interested in p°7r scattering, 
and only one of two decay channels of an I R — 1 resonance can contribute here, — > 
p ^ but not p ± 7T°. The second factor corresponds to the assumption that all three pion 
isospin states are equally populated in the gas. Similarly, for I R = the isospin factor is 
(1/3) x (1/3) = 1/9. We take into account the following resonances [9]: ai(1260), 7r(1300), 
a 2 (1320) and cj(1420). The nearest resonance under the threshold, a; (782), contributes a 
negligilbe amount due to its narrow width. For the spin factor we take F s — l, 1, 2, 1, 
respectively, for the aforementioned resonances. (These factors correspond to transverse 
p-mesons only.) The amplitude in the pion rest frame is obtained from Eq. (4) by 
multiplication by the rescaling factor k p /q cm = y/s/m n , where k p = \JE 2 — m 2 is the p 
momentum in the pion rest frame. 

For <7 pn we use the standard resonance formula. 
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. , 7T B R Yji 

r vim (v« - m ) + ry4 

According to Adler's theorem the pion scattering amplitude on any hadronic target van- 
ishes at zero pion energy in the target rest frame in the limit of massless pions. In the 
framework of an effective Lagrangian this can be achieved if the pion field enters by its 
derivative dip/dx p . We assume that in pir scattering through the a\ resonance dip/dx p is 
multiplied by the p-meson field strength tensor F pu and the a lu field. This results in the 
appearance of an additional factor in Ref pn and in a p7T in comparison to Eqs. (4) and (6) 



s — m 2 p — m% 
™ 2 ai -m 2 p -ml / 



(7) 



Here normalization at s — m 2 ai was performed. When s > m 2 ai this factor is replaced by 
unity. The analogous factors were also introduced for other resonance contributions. 
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At high energies we assume that the Regge approach is valid for 77? scattering and 
apply the vector dominance model (VDM) to relate pn and jn amplitudes. As is well 
known the Regge pole contributions to the forward scattering amplitude, normalized 
according to Eq. (2), have the form: 



/(«) = 



4ns 



, 1 



-ITTCti 



sin na,; 



-s ai ri 



(8) 



where ctj is the intercept of the i'th Regge pole trajectory, is its residue, and k is the 
projectile momentum in the target rest frame. As follows from Eqs. (2) and (8), 



(7(a) = £ r^" 1 , 



(9) 



Re/(s) 



4ns 



E 



1 + cos no.; 



-TiS 



Sin 7TQ!,- 



(10) 



For <r 77r only P (Pomeron) and P' Regge poles contribute [11,12]. The residues of the P 
and P' poles in 771 scattering were found by Boreskov, Kaidalov and Ponomarev (BKP) 
[12] using Regge pole factorization and data on jp, np and pp scattering. Taking BKP 
values of P and P' residues we have 



a WJ {s) 



7.48a 



ap—l 




a-pi — 1 



where a P = 1.0808, a P < = 0.5475, a = 1/137, s = 1 GeV 2 and a in Eq. (11) is given 
in millibarns. For P and P' intercepts we take Donnachie-Landshoff values [13]. Since 
in their fit of the data BKP assumed ap = 1, a pi = 1/2, the values of the residues in 
Eq. (11) are slightly changed in comparison with [12] in order to give the same value of 
cr 7r7 at s = 9 GeV 2 . From Eqs. (10) and (11) the real part of the forward 771 scattering 
amplitude can be found: 



k 

^&firy(,S)iriestbaine 7.48d; 



(\ ap — l / \ a P i— 1 



(12) 



where the momentum k is in GeV and Re/ is given in mb-GeV. 
In VDM cr p7T (s) is related to <r 77r (s) by [2] 



4na 



1 + -j) cr 77r (s), 



(13) 



where g 2 p /4n = 2.54, g 2 p j ' = 1/8, and the 99-meson contribution is neglected. A similar 
relation holds for Ref p7T . Unlike [2], we prefer here to use direct Regge formulae for Re/ 
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at high energies instead of inferring it from a by the dispersion relation since, in the latter 
approach]], the results are sensitive to the low energy domain, which is more uncertain. 

The results of the calculations of cr p7r and Re/ p7r as functions of p-meson energy in the 
pion rest frame are presented in Fig. 1. As may be seen from the figure the matching of 
low and high energy curves is satisfactory. 

For the amplitude Ref p ^ at laboratory energies of the p above 2 GeV we use the 
results of Ref. [|2| obtained with the dispersion relation, VDM and experimental data on 
<7 7 at. At lower energies we again use the resonance approximation 

The factor of 1/4 appears because we consider only transverse p-mesons. The isospin 
factors Fi are 1/3 and 2/3, respectively, for N and A resonances. We take 10 N and A 
resonances with significant branchings into pN and with masses above the pN threshold 
and below 2200 MeV as quoted in [9]. This set of baryonic resonances is close to the 
set used in [14]. The main difference in comparison with [14] is that the effective widths 
rpTV = r p V(q , cm/Q , ^ n ) 2i+1 were introduced only for the resonances close to the pN threshold 
(q^ is the value of the cm. momentum at the resonance). When q cm > q^ we put 
r^y- = T^ N . Besides these resonances, two others with masses below the pN threshold 
were accounted for: the A(1238) and the iV(1500). It was assumed that VDM is valid 
for the contribution of these resonances to the widths T pN and T^ N in the following 
form. Since both resonances are close to pN threshold, we can write for each of them 
Tpjv = QcmlpN and T^ N = k cm ^ lN) where q cm and k cm are the pN and momenta in the 
cm., respectively. Then we assume that 7 p at and 7 7 tv are related by the VDM formula 

7 7Ar = 47ra^r( 1 + ^ ]7 p7V . (15) 
9p\ gl) 

The value of 7 7 tv can be found from the values of o" 7 tv at the resonance peaks. The 
contribution of the A(1238) and of the iV(1500) to Kef P N are essential at low energies: 
they contribute about —1 to —0.5 fm at E p = 1 — 2 GeV in the nucleon rest frame. 

The results for <r pN and Ref P N in the rest frame of the nucleon, the curve obtained in 
Ref. 0] for high energies, and the matching curve are shown in Fig. 2. As can be seen the 
matching of low energy and high energy curves is good. 



■""We use this occasion to correct the misprint in the corresponding equation in Ref. [2] - in the 
equation (7) of [2] instead of the factor l/(27r) 2 it should be 1/27T 2 . In the calculations of [2] in fact the 
correct factor was used. 
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3 Determination of p-meson mass shift and width 
broadening based on the nucleon and pion distri- 
butions produced in heavy ion collisions. 

As mentioned above, in heavy ion collisions only nucleons and pions are considered as 
constituents of the medium. Therefore, in this case Eqs. (1) and (3) take the form 

Am(E) = [p N Ref pN (E) + p n Ref p7T (E)} , (16) 

AT(E) = - [p N a pN (E) + p w a p7T (E)} , (17) 

777- 

where p^ and p n are the nucleon and pion densities during the final stage of evolution of 
the hadronic matter produced in heavy ion collisions. 

We will restrict ourselves to consideration of central, head-on, collisions with small 
values of impact parameter when the number of participants - the nucleons, which undergo 
significant momentum transfer - is close to the total number of colliding nucleons. 

As shown by the experimental data, the nucleons and pions produced in heavy ion 
collisions cannot be considered as a gas in global thermal equilibrium even during the last 
stage of evolution of hadronic matter created in the collisions. In order to demonstate this 
let us discuss separately the cases of high energy, E ~ 100 GeV-A, and low energy, E ~ 
1 — 10 GeV-A, heavy ion collisions. In the case of high energy collisions the longitudinal 
and transverse momenta of nucleons and pions are very different. In the experiment on 
S + S collisions at 200 GeV- A [15] it was found that = 3.3 GeV, (ptn) = 0.61 GeV, 

and ~ 0.70 GeV, {pt-k) ~ 0.36 GeV. In other experiments on high energy heavy ion 
collisions - see [16,17] - the situation is qualitatively similar. This means that one can 
by no means speak about a thermal gas of final particles in global equilibrium, and their 
momentum distributions must be taken from experiment. This, however, still leaves open 
the possibility of local thermal equilibrium. 

The data for low energy heavy ion collisions also indicate that pions and nucleons 
cannot be described as gases in global thermal equilibrium. The angular distributions 
of pions produced in Ni + Ni collisions at E — 1 — 2 GeV A shows essential anisotropy 
[18]. If the pion angular distribution in the centre of mass system is approximated by 
1 + a cos 2 9 then, from the data, follows a 1.3. Unfortunately, there is not enough 
experimental information on nucleon angular and momentum distributions. We have 
checked the hypothesis of global thermal equilibrium by assuming that the probability 
of production of the given number of particles is proportional to the statistical weight of 
the final state ( Fermi- Pomeranchuk approach [19,20]). It is evident that this hypothesis 
is even more general than the hypothesis of global thermal equilibrium. In this approach 
the pion/nucleon ratio R n = n n /N in central collisions can be predicted in terms of the 
main ingredient of the method - the volume per particle at the final stage of evolution 
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and, of course, the initial energy. A calculation shows that the data [18] on the energy 
dependence of are well described by the statistical model, but in order to get the 
absolute values of R w in Ni + Ni as well as in Au + Au [21] collisions it is necessary to 
put the volume per nucleon very small, about 5 times smaller than in a normal nucleus, 
which is unacceptable. Therefore, the only way to perform the averaging over momentum 
distributions of pions and nucleons is by taking the latter from experimental data on 
heavy ion collisions. 

When calculating the p-meson mass shift and width broadening an averaging must 
be performed over the p-meson direction of flight relative to nucleons and pions. Such a 
calculation can be done only for real experimental conditions. For this reason we restrict 
ourselves to rough estimates. 

Consider first the case of high energies. As an example take the experiment [15] for 
central collisions where the ratio of pion to nucleon multiplicities was found to be R w = 5.3. 
Suppose that in this experiment the p-meson is produced with longitudinal and transverse 
momenta in the laboratory system ki = 3 GeV, kx = 0.5 GeV. We choose these values 
as typical for such an experiment. For these values of p momenta the formation time of 
the p-meson is close to the mean formation time of pions so a necessary condition of our 
approach is fulfilled. Since the mean momenta of nucleons and pions in the experiment 
[15] are known (they were presented above) it is possible, using the curves of Figs. 1 and 
2, to calculate the mean values of Ref P N, Ref p7T , a P N and <r p7T in pN and pn scattering. 
The results are, in lab frame: 

(Ref pN ) « -1.1 fm, (Re/ p7r ) « 0.03 fm, (18) 

(a pN ) « 45 mb , {a, m ) « 20 mb . (19) 

The small value of (Ref pn ) arises from a compensation of positive and negative contri- 
butions from low and high energy collisions, that is, from the scattering of the p-meson 
on pions moving in the same direction (comovers) or in the opposite one. Because of 
this compensation (Ref pn ) is badly determined, but since it is small this fact does not 
influence the final result. 

Using Eqs. (18) and (19) we can now find the mass shift and width broadening of the 
p-meson. For nucleon and pion densities we take 

N _ N 1 
PN ~V~Nv N + nv n ~ v N (l + R~^) ' ( } 

n w ^ Rw jr.-. \ 

Pw ~ V ~ Nv N + nv w ~ v N (l + R^) : [ ' 

where iV and n are the numbers of nucleons and pions at the last stage of evolution, 
R n = n/N, and V is the volume of system at this stage. It is assumed that at this 
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stage of evolution any participant - nucleon or pion - occupies the volume vn or 
respectively. We can write 

~ TTrj3 ' ^"ITi^' (22) 

where p° N = 1/fjv and (3 = v t /vn- For numerical estimates we take p% = 0.3 fm -3 , about 
two times standard nucleon density. This number is probably one of the most uncertain 
ingredients of our calculations. Substitution of Eqs. (18), (19) and (22) in (16) and (17), 
together with the experimental values R n = 5.3 and (3=1, gives 

Am p = 18 - 2 = 16 MeV , (23) 

Ar p « 150 + 400 = 550 MeV . (24) 

The first numbers above refer to the contributions from p — N and second from p — n 
scattering. Because the p-meson width broadening appears to be very large, a basic 
condition of our approach, Ar <C m p , is badly fulfilled. The applicability condition of the 
method, | Re/ |< d, is not well satisfied either since in this case d = 0.9 fm. For these 
reasons the values of Am p and Ar p may be considered only as estimates. 

The main conclusion to be drawn from Eqs. (23) and (24) is that for p-mesons pro- 
duced in high energy heavy ion collisions with the above chosen values of longitudinal and 
transverse momenta, the mass shift is small, but the width broadening is so large that 
one can hardly observe a p-peak in e + e~ or p + p~ mass distributions. Let us estimate 
how sensitive the results are to variations of and kx- It can easily be seen that the 
mass shift will always be small, say Am p < 50 MeV. If we put kx = instead of kx = 0.5 
GeV, this will only weakly influence the mean value of <j p n and decrease a pn by 20%. The 
latter results in a decrease of Ar p by 80 MeV, which is within the limits of accuracy of 
our estimates. The variation of k^ in the range 1 GeV to 10 GeV also results in variations 
of 10-20% in Ar p . 

As mentioned above, the main uncertainty in our approach comes from the assumed 
value of the nucleon density at the final stage of evolution: p° N = 0.3 fm -3 . If this 
density would be a factor of two smaller then Ar p ~ 250 MeV and the p-meson could be 
observed as a broad peak in the e + e~ or p + p~ mass spectrum. It should be mentioned, 
however, that the chosen value of (3 = v n /v]y = 1 is rather uncertain. If we assume that 
(3 = (tv/Vat) 3 , where r % and are pion and nucleon electromagnetic radii, 7v = 0.66 fm, 
tn = 0.81 fm, then (3 ~ 0.55. The choice of such (3 increases Ar by the factor of 1.6. 

In the course of p-meson propagation in the medium its decay width T(p — > tttt) may 
decrease?]. This effect can be estimated by substitution of an effective pion propagator in 
the medium [ k 2 — (m n — iT^/2) 2 ] _1 into the imaginary part of pion loop determining 

2 One of the authors - B.I. - is grateful to G.Brown for this remark. 
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p — > 7T7T decay. (Here r„- is the effective pion width in the medium arising from pion 
interaction with medium constituents.) The calculation performed in this way gives 



Even when T w ~ 500 MeV this correction is small. 

Our qualitative conclusion is that in central collisions of heavy nuclei at high energies, 
E ~ 100 GeV-A, where a large number of pions per participating nucleon is produced, 
the p-peak will be observed in e + e~ or p + p~ mass distributions only as a very broad 
enhancement, or even no enhancement at all. Inspite of the assumptions we made, in- 
cluding noninteracting nucleon and pion matter at the final stage of evolution and the 
specific numerical value of the nucleon density, we believe that this qualitative conclusion 
is still valid. This conclusion is in qualitative agreement with the measurement of e + e~ 
pair production in heavy ion collisions [22] where no p-peak was found and only a smooth 
e + e~ mass spectrum from to 1 GeV was observed. If, however, such a peak would be 
observed in future experiments it would indicate that the hadronic (nucleon and pion) 
density at the final stage of evolution, where the p-meson is formed, is very low, even 
lower than normal nuclear density. 

Recently preliminary data in Pb — Au collisions at 160 GeV- A have been presented 
[23] where, in studying the e + e~ mass spectrum, it was found that the p-peak is absent 
at kr(e + e~) < 400 MeV, but reappears as a broad enhancement at A>r(e + e~) > 400 MeV. 
We do not see the possibility for such a phenomenon in the framework of our approach for 
central heavy ion collisions. Moreover, we believe that for central collisions the absence 
of a p-peak at low kx and its reappearance at higher kx will be hard to explain in any 
reasonable model. The only explanation we see for this effect is that in this experiment 
peripheral p-meson production plays an essential role. Then p-mesons with higher kr 
have a larger probability to escape the collision region and decay as free ones. 

Let us turn now to the case of lower energy heavy ion collisions, E ~ a few GeV-A. 
Consider, as an example, heavy ion collisions at E^ in — 3 GeV- A and production of p- 
mesons of energy i?* ot = 1.2 GeV in the forward direction. (This particular value of 
the p-meson energy was chosen because our approach works better at higher E p , and 
p-mesons of this energy can be kinematically produced at such a beam energy). The 
number of pions produced can be estimated by extrapolation of the data [18] on Ni + Ni 
collisions. This data shows, with good accuracy, that R n is linear in \fs/2 — m. We find 
that R n = 0.48. As follows from analysis of the data [18] at E^ n = 1.93 GeV-A, the 
average energies of produced pions are rather small: E w ~ 200 — 300 MeV. At such low 
energies it is reasonable to suppose that for pions (p L ) = (p±) 0.2 GeV. Assuming that 
the mean perpendicular momentum of nucleon participants is the same as at high energy 
- (ptn) — 0.61 GeV [15] (this assumption does not much influence the final results) we 
can construct the momentum distributions of nucleons. Then we are in a position to 




(25) 
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calculate the mean values of Ref P N, Re/ pvr , <j p n and a p7T . The results are: 

(Ref pN ) = -0.54 fm , (Re/ p7r ) = 0.30 fm , (26) 

(o~pn) = 45 mb , (cV-) = 13 mb . (27) 

For the p-meson mass shift and width broadening we have, with the same value of p% as 
above and (5 — 1: 

Am p = 37- 10 = 27 MeV , (28) 

A T p = 245 + 35 = 280 MeV . (29) 

The first numbers above refer to pN scattering, the second ones to pn. The conclusion is 
that in low energy heavy ion collisions a p-peak may be observed in e + e~ or p + p~ mass 
distributions as a broad enhancement approximately at the position of p-mass. 
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Figure Captions 



Figure 1: Cross section (a) and real part of the forward scattering amplitude (b) for p- 
mesons scattering on pions as functions of the total p-meson energy in the pion rest frame. 
The curves at low energy are the result of the resonance approximation. The curves at 
high energy are the result of the Regge parametrization. These curves are matched at 
intermediate energies. 

Figure 2: Same as for Fig. 1 but for p- mesons scattering on nucleons. The curves at high 
energy are from Ref. [2]. 
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